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Abstract 

Wc extend the work o f lDelong and Imkellerl (j2010allbh concerning Backward stochastic differential 



equations with time delayed generators (delay BSDE). Wc give moment and a priori estimates in 
general L p -spaces and provide sufficient conditions for the solution of a delay BSDE to exist in 
LP . We introduce decoupled systems of SDE and delay BSDE (delay FBSDE) and give sufficient 
conditions for their variational differentiability. We connect these variational derivatives to the 
Malliavin derivatives of delay FBSDE via the usual representation formulas. We conclude with 
several path regularity results, in particular we extend the classic L 2 -path regularity to delay 
FBSDE. 
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Introduction 

The theory of nonlinear backward stochastic differential equations (BSDEs) wa s introduced by 
Pardoux and Peng ( 1990l ) with its main mo tivations being mathem atical finance (see El Karoui et al. 
( 19971 )) and stochastic control theory (see Yong and Zhou ( 19991 )). In the last twenty years much 
effort has been given to this type of equations and nowadays many classes of BSDEs and results 
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on them are available. Due to tractability, common results are achieved within a Markovian frame- 
work. Under certain conditions the BSDE's solution exhibits a Markov structure and hence can be 
interpreted as an instantaneous transformation of the underlying Markov process that spans the 
stochastic basis of the underlying probability space. This in turn yields access to the theory of 
partial differential equations via the non-lin ear Feynman-Kac formula. 

Moving away from the Markovian setting, iDelong and Imkeller fcoiOal lbh introduce a new class 



of BSDE labeled backward stochastic differential equations with time delayed generators (delay 
BSDEs). The dynamics of these BSDEs are governed by 



Yt = £+ f f(s,Y(s),Z(s))ds- [ T Z s dW s , t£[0,T], 
Jt Jt 

where the generator / at time s € [0, T] is allowed to depend on the past values of the solution (Y, Z) 
over the time interval [0, s] and £ is a measurable random variable. In these two works the authors 
answered thoroughly several fundamental questions: existence and uniqueness of a square integrable 
solution, comparison principles, existence of a measure solution, BMO martingale properties for the 
control component Z of the solution, Malliavin differentiability for delay BSDEs driven by a Wiener 
process and a generalized Poisson martingale. To the best of our knowledge the only existe nce and 



uniqu eness results for this class of BSDEs follow from those two works. As pointed out by IDelong 



(2010). delay BSDEs appear naturally in finance and insurance related problems of pricing and 



hedging of contracts. In the same work the author analyses a vast scope of contracts to which this 
class of BSDEs can be applied to. 

Paying consideration to and seeking reference from the state of the art of BSDEs with non-time 
delayed generators, the next step concerning delay BSDEs is to obtain a feasible numerical scheme. 
Here, the main obstacle is the presence of the control process Z in the generator. This process is 
usually obtained via the predictable representation property of the underlying stochastic basis, and 
initially all one knows about Z is that it is a square integrable process. To steer in the direction 
of a numerical scheme a deeper analysis on the fine properties of the solution of such equations is 



requir ed. As for numerics fo r Lips chitz continuous BSDEs (see for example iBouchard and Touzi 



t()04 ) or lBender and Zhang] fcOQSl )) one is usually forced to gather several results concerning the 



path regularity properties of the solution process before being able to give proper convergence 
results. Such path properties include not only sample path continuity but also estimations on the 
time increments of the components of the solution by the size of the time increment. For the purpose 
of establishing such path properties we first need to prove several auxiliary results. 

Our agenda consists of refin ing and extending the existence and uniqueness results obtained in 
Delong and Imkeller! (|2010al lbh and then steer into the direction of the smo othness properties of 



the sol ution of delay BSDEs. We start by improving the original results of IDelong and Imkeller 
(|2010ah concerni ng their a priori estimates by reformulating them in a more standard fashion. In 
Lemma 2.1 from IDelong and Imkellerl ( 2010al ). the a priori estimates expresses the difference (in 



norm) of the solution of two delay BSDE as the difference of the respective terminal conditions and 
generators. These a priori estimates fall short of the usual a priori estimates one expects to see due 
to the presence of the solutions of both delay BSDE on the right hand side of the estimate. We 
establish a priori estimates in the classical form where the right hand side of the estimate contains 
the difference of generators evaluated at their zero spatial state and hence i s independent of the 



BSDE solutions. Within the topic of a priori estimates we extend the results of lDelong and Imkeller 
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(l2010ah in another direction. We show that given extra integrability of the terminal condition and 
the generator, the solution will inherit this integrability. This allows us to state moment and a 
priori estimates in general L p -spaces and no t solely in L? . The proof of these estimates relies on 
techniques from iDelong and Imkel ler (2010a) and on computations carried out for non-time delayed 



BSDEs in the spirit of IWang et al.l ([20071). The usual techniques to obtain higher order moment 
estimates fail in the setting of delay BSDEs, the reason for this will be seen in (jlip below. A rough 
explanation would be that for the usual (non-delay) BSDE setting the dynamics of Yj is given by 
sums of Lebesgue and Ito integrals over the interval [t,T] but for delay BSDEs the dynamics of Yt 
depends also on a integral over the whole interval [0, T] which doesn't allow the usual techniques 
to be used. The general estimates we obtain pave the way to a result of existence and uniqueness 
of solutions to delay BS DE with Lipschitz continuou s generators in general LP spaces for p > 2. 
Inevitably, in analogy to Delong and Imkeller ( 2010al fbl) a compatibility condition on the Lipschitz 
constant and terminal time is required to obtain existence of solutions (see our Theorem 12. 14j) . 

A customary field of application of BSDEs consists in coupling them with SDEs, giving rise (in our 
case) to systems of delay forward-backward SDEs (delay FBSDEs). We show that when coupling 
a delay BSDE with a forward diffusion and assuming appropriate regularity conditions, we obtain 
smoothness properties of the solution in terms of the involved parameters, in particular with respect 
to the init ial condition of the forward d iffusion. Combining this with the Malliavin differentiability 
proved in Delong and Imkeller ( 2010bl ) enables us to derive the usual representation formulas for 
FBSDE which display the relationship between the Malliavin derivatives of the solution process and 
their variational (classical) derivatives. It is somewhat surprising that such a relationship still holds 
since it is usually a consequence of the BSDE's Markov property which clearly fails to materialize 
in the context of delay FBSDE. 

With this collection of results we are finally able to address the p ath regularity issue of delay BSDE. 
Using the techniques employed in Imkeller and Dos Reisl ( 2010al lbh . we establish path continuity for 
the components of the solution of delay FBSDE and we give a result that bounds the norm of the 
increments in time of Y and Z by the size of the time increment. We expect that these results 
will open the door to the derivation of concrete numerical schemes and their convergence rate and 
intend to tackle these problems in our future research. 

The paper is organized as follows: in Section [1] we fix notations and elaborate on the type of 
time-delaye d BSDEs that we co n sider. In Section [2] we refine and extend the a priori estimates 
obtained in lDelong and Imkellerl (|2010al ) and then use them to establish existence and uniqueness 
of solutions in general LP spaces. In Section [3] we introduce the delay FBSDE framework and use 
results from the previous sections to obtain the differentiability of the solution process with respect 
to the initial state of a forward diffusion. The representation formulas and the path regularity 
results are presented in Section [H 



1 Preliminaries 

Let (fi, J-, P) be a probability space equipped with a standard cf-dimensional Brownian motion W. 
For a fixed real number T > we consider the filtration F := (J-t)t>o generated by W and augmented 
by all P-null sets. The filtered probability space (fl, J-, F, P) satisfies the usual conditions. Depending 
on whether we work on M. d or M. mxd , the Euclidean norm respectively the Hilbert-Schmidt operator 
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norm is denoted by | • |. Furthermore, V denotes the canonical gradient differential operator and 
for a function h{x,y) : W 71 x M. d — > W 1 , we write V x /i or V y h for the derivatives with respect to x 
and y. We work with the following topological vector spaces: 

• For p > 2, let L p (IR m ) be the space of J^-measurable random variables £ : f2 — > M m normed 
byU\\ LP :=E[|^] 1/P . 

• For p > and p > 1, H p JR m><d ) denotes the space of all predictable process ip with values in 



ornxd 



such that the norm ||y||-^p := E ^ e^ s \ip s \ 2 ds S j 



sp/2 



l/p 

< oo. 



For p > and p > 2, S p JM. mxd ) denotes the space of all predictable processes r\ with values 



in M mx " such that the norm H^H^p := E 



sup 0<t <T e pt \rj t \ 2 



p/2- 



i/p 

< oo. 



We omit referencing the range space if no ambiguity arises. It is fairly easy to see that for any 
P, P > the norms on H^, H 1 ^ and 5^, are equivalent. 

Some notation 

We introduce a notational convention which will be used throughout the text: for an arbitrarily 
given integrable function / : [0,T] — > W n , trivially extended to [— T, 0) via /(t)lr_r i0 )(*) = 0) an d 
a given deterministic finite measure a supported on [-T, 0) which is not necessarily atomless, we 
denote for t £ [0, T] and any p > 2 

(f ■ a)(t) := f° f(t + v)a(dv) and (f ■ a)(t) := f° \f(t + v)\ p a(dv). 
J-T J-T 

Similarly, for a given process {}Pt)te\o,T\i extended to [— T, 0) by imposing ip t = on [— T, 0), we 
denote 

(<p-a)(t) := / Vt +V a(dv), t € [0, T], (1) 



and 

f0 



(ifP-a)(t):= [ \<p t+v \ p a(dv), te[0,T],p>2. (2) 

J-T 

We now give a lemma concerning the change of integration order for ([I]) and ([2]) , which will become 
useful in the sequel. 

Lemma 1.1. Let ip be a process and a a non-random finite measure supported on [— T, 0). Then 
we have the following change of integration order: for every k > 1 

T r-T 

(ip k ■ a)(s)ds = / a([r-T, (r - t) A0))\ip r \ k dr, Vt G [0, T], P — a.s. 
Jo 

Moreover, if we have for p > 1 £/ia£ y € %q, i/ien we a/so /iaue t/iai 

||(^-a)||^<M p |M|^, 

w/iere M p = (e /3T ) p / 2 (a([-T, 0))) p . 
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Proof. Let t in [0,T] and k € [1, +00). We have that 



T 



(tp ■ a)(s)ds 



T r 



t J-T 
[-T+v 



\ip s+v \ a(dv)ds 



rT 
-T Jt 



(p s+v \ h ds a(dv) 



T r(l — t)AO 







\(p r \ dr a(dv) 

T J(t+v)\/0 
T 

a([r — T, (r -t) A 0))\ip r \ k dr. 



J(r-T) 



\ip r \ k a(dv) dr 



The second claim follows by applying Jensen's inequality and changing the integration order as done 
above, i.e. for any (3 > and p > 1 we have 



E 



T \ p/2 

e^ s \(<p ■ a)(s)\ 2 ds 



< (e^ T a([-T,0))) p/2 E 

\P/2 



(\^-a)(s)ds) 



\P/2 



< M p E 



which concludes the proof. 



□ 



2 General results on BSDE with time delayed generators 

In this section we give a brief overview of BSDEs with time delayed generators and discuss the 
setting they are studied under. We then establish convenient a priori estimates on the difference of 
two solutions to such equations which will play a central role in proving existence and uniqueness 
of solutions in the more general W-spaces. 

2.1 BSDEs with time delayed generators 

Let us start with a recap on BSDE with time delayed generators. Throughout the text, we assume 
(HO) a y ,a z are two non-random, finitely valued measures supported on [— T, 0) 
We also define 

a:=a y ([-T,0))y a z ([-T,0)). (3) 
Given p > 2, we assume that the following holds: 

(HI) £ is an J^-measurable random variable which belongs to L p (M m ); 

(H2) the generator / : O x [0,T] x M. m x R mxd — > M, m is measurable, F-adapted and satisfies the 
following Lipschitz like condition: there exists a constant K > such that 

\f(t, y, z) - f(t, y', z')\ 2 <K(\y- yf + \z - zf) 

holds for dP ® dt-almost all (co, t) € Q x [0, T] and for every (y, z), (yV) £l m x R mxd ; 
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(H3) E[(/ T |/( S ,0,0)| 2 d S ) p/2 ] <oo; 
(H4) f(t,-,-) = if t <0. 

Following the notation from equation ([I]), we write 

(Y-a y )(t) = / Y t+V a y {dv) and (Z ■ a z )(t) = / Z t+V a z (dv), 0<t<T, 

J-T J-T 

for some processes (Yt)te[o,ri an d (Zt)te.[o,T] satisfying appropriate integrability conditions. Assump- 
tion (H2) and Jensen's inequality then imply 



(H2>) \f{t, (Y ■ a y )(t), (Z ■ a z )(t)) - f(t, (Y' ■ a y )(t), (Z' ■ a z ){t)) \ 2 
< K{\((Y - Y') ■ a y )(t)\ 2 + \((Z - Z') ■ a z )(t)\ 2 } 
<L{{{Y-Y'f.a y ){t) + {{Z-Z'f 



a.2 
«*)(*)}> 

where L := Ka with the real number a given by (|3|). The focus of our study are BSDE with time 
delayed generators which are of the type 



Yt = £+ [ f{s,T(s))ds- [ Z s dW s , 0<t<T, 
Jt Jt 



where V abbreviates for t € [0, T] 

r0 



T(t) :- 



Y t+V a y (dv), J Z t+V a z (dv)j = (jY ■ a y )(t), (Z ■ a z )(t)j . 



(4) 

(5) 
9 if 



Definition 2.1 (Solution of a Delay BSDE). We say (Y,Z) is a solution to the delay BSDE 
(Y,Z) belongs to the space Sq x Hq and satisfies (jU). 

Using a fixed point argument, Delong and Imkeller ( 2010al ) have shown that a BSDE of the type 
(j3])-([5]) admits a unique solution if the parameters of the equation (jl]) are sufficiently small, i.e. 
if the Lipschitz constant K > or the terminal time T > satisfy a smallness condition. The 
following L 2 -existenc e and u niqueness result is a straightforward modification of Theorem 2.1 from 



Delong and Imkellerl (|2010al ). 



Theorem 2.2. Let p = 2 and assume that (H0)-(H4) are satisfied. For a defined as in ([3]), assume 
that the non-negative constants T, L = Ka, (3 are such that 

1 f° 

(8T + -)L e- /3lt p(dti)max{l,T} < 1, for p G {a y ,a z }. 
P J-T 



Then the delay BSDE ([I])-© has a unique solution (Y, Z) G 5|( 

Remark 2.3. In \Delona and lmkelle\ koiod ). this result is proved for the one- dimensional case 
d = m = 1. It is clear that by the nature of the fixed point argument, the proof is insensitive to the 
dimension of the equation. 

Remark 2.4. Given that a compatibility condition is necessary in order to establish existence and 
uniqueness of solutions and moreover that we will be giving an extended version of it, all the proofs 
in this section are given with extra detail in order to better control the constants involved in each 
result. 
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2.2 Moment and a priori estimates 



In Lemma 2.1 from iDelong and Imkellerl (|2010al ) the authors provide a priori estimates for the 
time delayed BSDE @ which estimates the norms of the difference between the solution of two 
BSDE in terms of the terminal condition and the difference of the generators applied to the solution 
processes. More specifically, for i G {1, 2} let (Y l , Z l ) be the solution of a BSDE with dynamics (jH) 
with terminal condition £ l and driver f l satisfying (H1)-(H4), then it holds that 



|Y J 



Y 2 \\ 2 n2 + 



iz 1 



Z 2 \\ 2 H 2 < C 2 



E[ 



2|2 " 



+ e?'\f\8, (Y 1 ■ a)(s), (Z 1 • a)(s)) - f 2 (s, (Y 2 ■ «)(«), (Z 2 • a)(s))\ 2 ds] }, (6) 

where the autho rs assume that a is some de terministic measure on [— T, 0) with mass one. Thus 



Lemma 2.1 from lDelong and Imkellerl (|2010al ) establishes the a priori estimate Q whose right hand 
side de pends again on the solution of both delay BSDE. In the context of IDelong and Imkeller 
l|2010aJ ) such a result suffices to establish existence and uniqueness of solutions in iS 2 , x H 2 but 
the situation becomes more intricate when the same issues are considered on S p x Ti^ for p > 2. 
More precisely, we are not able to obtain an estimate similar to © when p > 2. In addition, the 
study of differentiability of the solution (for both p = 2 and p > 2), made in Section [31 requires a 
priori estimates where the right hand side of the estimate depends only on the problem's data: the 
differences between the terminal conditions and a quantity of the form 82 f s '■= f 1 {s, (Y 2 ■ a)(s), (Z 2 ■ 
a)(s)) — f 2 (s, (Y 2 -a)(s), [Z 2 -a)(s)). For a clear view of the required estimates, compare for instance 
(H with flSb. 



Moment estimates - part I 



x U P p(R mxd 



As a starting observation, we have that if ([!]) admits a solution (Y, Z) in H 
then we also have that Y G S^(R m ). 

Lemma 2.5. Let (3 > 0, p > 2 and assume that (HO)-(H^) hold. If the delay BSDE @ admits a 
solution (Y,Z) G ^(M m ) x 7^(R mxd ) f/jerc we /iaue a/so t/iat Y G 5^(M m ). 

Proof. Throughout let t G [0, T] and p > 2. Since all /3-norms are equivalent, it suffices to show the 
result for /3 = 0. We drop the /3-subscripts in the following. The pair (Y, satisfies 

y t = e+ / T /(s,(y-a y )(s),(^-a^)(s))d s - [ T Z s dW s , 

Jt Jt 

and in turn we have 

sup |y t |<|e|+ [ T \f{s,(Y-a y )(s),(Z-a z )(s))\ds+ sup | F Z,dW a \. 

0<t<T JO 0<t<T Jt 

Combining the fact of Z G W with the inequalities by Young, Doob and Burkholder-Davis-Gundy 
(BDG), we obtain 



E 



sup 

0<t<T ' Jt 



T 



ZAWs 



p/2- 



< 2 p/2 E 

< 2 P E 



z s dw s r+ sup / z s dw s 

0<t<T Jo 



2 y/2 



sup 

0<t<T JO 



Z,dW,\ p 



< 2 p C p \\Z\\P < 00. 
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Next observe that by the Lipschitz property of the generator / (notice that (H2) implies (H2')), it 
follows that 



T 



9 \p/2 

f(s,(Y-a y )(s),(Z-a z )(s))\ 2 ds) 



2 xp/2 

da) 



<2 p/2 (j o \f(s,0,0)\ 2 ds + f(s,(Y-a y )(s),(Z-a z )(s))-f(s,0,0) 

< 2P/ 2 2^ 2 - 1 |(^ T |/( S ,0,0)| 2 d S ) P/2 + (l£ ((|Y| 2 • a y )(s) + (\Z\ 2 ■ a z )(s))ds) 
The second term in the bracket can be further estimated by 

(l£ ((|Y| 2 • a y )(s) + (\Z\ 2 • a z )(s))ds *' 



\P/2 



< 2 p/2-l L p/2 



< 2 p/2_1 L p/2 a p/2 



T 



l*T-a- 



T 



,)(s)ds) 



\P/2 



+ 







\Z\ 2 ■ a z )(s)d s y /2 



,9 W 2 / f , ,9 \ P / 2 

y s | 2 ds) + ( / |z,| 2 ds) 



where the last line follows from Lemma ll. 11 This estimate together with (H3) yields 



E 



T 



|/( S ,(y-a y )( S ),(Z-^)( S ))| 2 d S 



p/2 



< oo. 



Using hypothesis (HI), i.e. that £ is in € L p , we can conclude that Y S S p must hold. 



□ 



A priori estimates 

Let us define the weighted variant a of a as the maximum of the weighted measures a y and a z on 
i-T,0) by 



rO 

e~ l3s a y {ds)y e-P s a z (ds), f3 > 0. 

-T 7-T 



(7) 



Remark 2.6. M^e emphasize that a depends on j3. To keep the notation to a minimum we simply 
write a instead of making the dependence explicit. 

The next results establishes canonical a priori estimates (in the sense that the right hand side of 
the estimate only depends on the problem's data) for the solutions of two time-delayed BSDEs as 
given by ([!]). We distinguish between the cases p = 2 and p > 2, and we start with the case p = 2. 

Proposition 2.7 (A priori estimates for p = 2). Let p = 2. Consider i E {1, 2} and let (Y l , Z l ) £ 
Sq x H 2 be the solution of the delay BSDE @ with terminal condition £ l and generator f l satisfying 
(H0)-(H4). Denote by K > the Lipschitz constant of f 1 as given in (H2 J ) and set 5Y = Y 1 — Y 2 , 
5Z = Z l — Z 2 . If either T or K or a are small enough then there exist two constants /3,7 > 
satisfying 

olL c\L 
D\ := j3 — 7 > and D2 := 1 > (with L = Ka and a as in flSJ)), (8) 

7 7 
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and a constant C 2 = C 2 (/3, 7, a, L, T) > depending on /3,^,a, L,T such that: for i G {1,2}, 

cT 



(J\Z { ) eSjxTil and 



\\SY\\ 2 sj + \\SYf n} + \\5Zf n} < C 2 {e^ t \6Y t 



+ E 







\hfs\ 2 &s\); 



(9) 



where 6 2 f t := f 1 (t, (Y 2 • a y )(t), (Z 2 • a y )(t)) - f 2 {t, (Y 2 • a y )(t), (Z 2 ■ a y )(t)) for t G [0, T]. 



Proof. Let 7, if, T, a be such that the relations in (JSj> are satisfied (i.e. D\ > and D 2 > 0). 
Throughout let t G [0, T], i G {1, 2} and define T l as in ([5]) for the pair (Y l , Z % ). An application of 
Ito's formula to the semimartingale e^*|5Yf| 2 for /3 > yields 

f T f T 

e? t \5Y t \ 2 + J f3e ps \5Y s \ 2 ds + J e? s \5Z s \ 2 ds 

= e P T \5Y T \ 2 + / T 2e^(«5y s ,/ 1 ( S ,r 1 ( S ))-/ 2 ( S ,r 2 ( S ))>d S - [ T 2e^(5Y s ,5Z s dW s ) 
Jt Jt 

< e? T \5Y T \ 2 + j\e? s \5Y s \ 2 ds + ^ — (\f(s, r\s)) - f(s, r 2 ( S ))| 2 )d, 

+ 2 f eP a (SY Si 6 2 f.)ds- [ 2e? s {6Y s , 8Z s dW s ), 
Jt Jt 

where the last inequality results from Young's inequality for 7. Reorganizing and taking condition 
(H2') for the generator f 1 into account, we get 



e^\5Y t \ 2 + [ T (f3- 7 y s \5Y s \ 2 ds+ f e^ s \5Z s \ 2 ds 
Jt Jt 



< e? T \5Y T \ 2 + / T e —L\{\5Y\ 2 ■ a y )(s) + (\5Z\ 2 ■ a z ){s) 
Jt 7 L 

r-T fT 



d.S 



+ 2/ e Pa (5Y a ,5 2 f 8 )ds- [ 2e" s (6Y s ,5Z s )dW s . 
Jt Jt 



By a change of integration order argument similar to that in the proof of Lemma 11.11 we obtain for 

j G {y, z} and <fP> = 5Y, (f> z = 5Z 



e^d^f ■a j )(s)ds 



T r 



T 

r-T+v 



^ {s+v) e^ v t {s+v > 0} \4 + v\ 2 ^(dv)ds 




< 



T J(t+v)V0 



3^e-^l {r > 0} |^| 2 dr aj (dv) 



T p(r~t)A0 




Jr-T 



e-^a. 



T 



j(dv))dr < / SeH<#l dr, 



(10) 
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with a given by (J7J). Continuing the inequality from above we get 



e^\5Y t \ 2 +[ (p- 1 ) e P s \6Y s \ 2 ds+ [ e? s \5Z s \ 2 ds < e? T \5Y T \ 2 + 2 / e^(5Y s ,5 2 f s )ds 
Jt Jt Jt 

+ [ ^ S (\5Y S \ 2 + \5Z s \ 2 y s - j\e? s (5Y s ,5Z s dW s ). 



(11) 



Taking the expectations for t = yields 



7 



T 



(0-7-— )E / e^\SY s \ 2 ds +(!-— )E / e^ s \5Z s \ 2 ds 



o 



7 



< E 

< E 

< E 



e^ T \5Y T \ 2 
e^ T \5Y T \ 2 



+ 2E 
+ 2E 



T 



e l3s (6Y s ,5 2 f s )ds 



e^ T \5Y T \ 2 + 7 'E sup e pt \8Y t f 
- 1 L 0<t<T 



sup e^\SY t \ / e * s |<$ 2 / fl |ds 

0<t<T Jo 

T 



+ ^tE 

7 



e2 



'|<52/ a |da) ; 



where we have used Young's inequality with some 7' > to be specified later. From the last 
expression and since D\ , D 2 > (see ([8]) ) we deduce that 

\\5Y\\ 2 n2 +\\5Z\\ 2 n2 <cU\e? T \5Y T \ 2 } +7'\\5Y\\h + A E [( T e^|5 2 /,|d S ) 2 ), (12) 

where C > is a constant depending /3,7, a,L and T. In order to obtain the iSj-estimate for <5Y 
we observe that we have 

5Y t <5Y T + [ \f l (s,T\s)) - f 1 (s,T 2 (s))\ds+ \6 2 f s \ds - F 5Z s dW s . 
Jt Jt Jt 

Multiplying by the monotone increasing function el* and taking the conditional expectation with 
respect to Ft we get 



e^SY t < E 



< E 



E 



e §t\5Y T \ + el* ^ \f l (s, T l (s)) - />, T 2 (s)) \ds + ef * ^ \5 2 f s \ds\F t 
e^ T \5Y T \ + j\^\f l (s,T\s)) -f l (s,T 2 (s))\ds 

- f e^ a \f 1 {s,Y\s)) - / 1 (s,r 2 ( S ))|d S + [ T e% s \5 2 f s \ds+ f e^ s \8 2 f s \ds\F t 
Jo Jt Jo 

e% T \5Y T \+ [ T e^ s \f l (s,T l (s)) - f 1 { S ,T 2 (s))\d S + C e^\S 2 f s \ds\T t 
Jo Jo 
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Using Doob's inequality, we obtain 



\SY\\% 



< 4 E 



< 12 E 



E 



e% T \6Y T \+ e? s \f 1 {s,T 1 ( S ))-f l (s,T 2 ( S ))\ds + e? s \6 2 f s \ds \ T T 



e^\5Y T \ 2 + T e^\f 1 (s,r\s))-f 1 (s,T\s))\ 2 ds + 



I, 
e2- 



\hfs\ds)' 



where the last line follows by Jensen's inequality. Since f 1 satisfies (H2'), an application of Lemma 
O yields 



||*y||^<12^E 



>/ T \5Y T \ 2 } +aTL(\\5Y\\ 2 nl + \\5Z\\ 2 nl )+E[{ f e% s \6 2 f s \ds) 2 ] }. 



Hence, plugging into (|12p we find 

(1 - 12C-/'aTL)E\ sup e pt \5Y t \ 2 



-0<t<T 

< 12{(l + CdTL)E 



e pi \5Y T \ 2 + 1 + Cy ctTLjE 



e^\S 2 f s \ds) 2 ]y 



Choosing 7' small enough such that (1 — 12C^f'aTL) > is satisfied we conclude that estimate ([9]) 
holds for a constant C 2 = C 2 (f3, 7, a, L, T). □ 

Remark 2.8. Note that in the previous result we have three degrees of freedom: the Lipschitz 
constant of the driver K, the time horizon T and the duration of the time delay given by a. 

The proof for the case p > 2 is more involved and uses techniques from the proof of Proposition 
12.71 The main reason for the proof to be more involved can be seen in (jlip . Usually the dynamics 
of Yf is described by integrals over the interval [t, T] but for delay BSDEs we see from (jlip that the 
dynamics of Y± depends also on a i ntegral over the whole interval [0, T\. We also remark that the 
techniques of Delong and Imkeller ( 2010al ) cannot be extended in LP (for n > 2 ), see for instance 
estimate (2.3) present in the proof of Lemma 2.1 in Delong and Imkeller! ( 2010a ). 

The next proposition gives a result that will be central in establishing existence and uniqueness of 
L p -solutions to delay BSDEs as well as in proving the differentiability results of Section [3j 

Proposition 2.9 (A priori estimates for p > 2). Let p > 2. Consider i € {1,2} and denote by 
(Y l ,Z l ) € Sq x T~Lq a solution of the delay BSDE @ with terminal condition and generator f % 
satisfying (H0)-(H4). Denote by K > the Lipschitz constant of f 1 in (H2') and set 5Y = Y 1 — Y 2 , 
5Z = Z l — Z 2 . Lf either T or K or a are small enough (for L = Ka, a as in ^ and a as in ([7]) ) 
then there exists /3,7 > satisfying ([8]) (i.e. D\, D 2 > 0) and 



1 



^ p -'d 2 p/2 \ 



P \p/2 



aL 



' ^ r\i r\)T ' 



p/2 -p/2 



f'V" 2 > (13) 



"p — 2' 7 — aL' 7 ' p — 2 

where m G N denotes the dimension of the 5Y process and the constant d p / 2 is given by 



d 



p/2 - z 



m 



p/2+1/ P ^P 2 /2fp{p-D\P/ 2 



(14) 
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In addition, (Y l , Z l ) E x (i € {1, 2}) and there exists a constant C p = C p ((3, 7, a, L, T, m) > 



explicitly given in (|26j) such that 

\\SYf sP + \\SY\\r + \\6Z\\ P nP < C p \e \(e^ T \5Y T \ 2 ) p/2 



+ E 



J e^\5 2 f s \ds) p ]}, (15) 



with S 2 f t = f l {t,Y 2 (t),Z 2 (t)) - f 2 (t,Y 2 (t),Z 2 (t)), forte [0,T\. 
Remark 2.10. A closer analysis on the constants D\, D 2 and D3 shows: 

lim (D U D 2 ,D 3 ) > (0,0,0). 

This means that with either a small T or a small K or a small a the conditions of the previous 
result can be verified. 

Proof of Proposition [OJ Throughout let t G [0, T], i 6 {1, 2} and from define D\ := (3 - 7 - 

and D 2 := 1 — We emphasize that a as defined in (J7|) depends on (3. Recall fill j) from the proof 
of Proposition 12.71 



|<5Y t | 2 + / T (/3- 7 )e^|<5y s | 2 d S + /V s |£Z s | 2 d S <^ T [SY T | 2 + 2 C e^ s {SY s , 8 2 f s )ds 
Jt Jt Jt 



aL 



(16) 



By assumption (3, 7, T, K, a are such that (JH]) holds and hence we have that D± > and D 2 > 0. 
We carry out the proof in several steps. 

Step 1: We claim that 



E 



V|^| 2 ds V 







< D, 



E 



" \$ y t\ ) 



2^p/2 



+ 2^ 2 d 2 p/2 D^ 2 \\5Y\\ p Pg 



+ 2 3,/2-i E [|^ e ^<5y s ,5 2 / s >d s | p/2 ]}, 



(17) 



where d p / 2 > is a given constant appearing in the BDG inequality which only depends on p > 2 
and the dimension. Estimate (JT7J) can be deduced as follows: putting t = in (|16p and noticing 
that by (|8|) the constants Z?i and D 2 are positive we get 



aL-, 



7 Jo 



aL 



7 Jo 



(l-— ) / e Ps\5Z s \ 2 ds<(P- 1 -—) e ? s \5Y s \ 2 ds +{!-—) e P s \5Z s \ 2 ds 



7 Jo 

T 



<e? T \5Y T \ 2 + 2 [ e^(5Y s ,5 2 f s )ds-2 [ e^ s (5Y s ,5Z s dW s ). 
Jo Jo 

Now raising both sides to the power p/2 > 1, making use of the fact that for a, b, c € M 

| a + 2b - 2c\ p/2 < 2 p l 2 ~ x (\a\ p ' 2 + \2b - 2cf / 2 ) < 2^' 2 " 1 (\a\ p l 2 + 2 p ' 2 ~ l (\2b\ p ' 2 + \2c 
= 2 p/2 ~ 1 \a\ p/2 + 2 3p/2 ~ 2 \b\ p/2 + 2 3p/2 ~ 2 \c\ p/2 
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and taking expectations, we get 



1 - ^V 2 E 
7 } 



V|<5Z s | 2 d S V 



< 2 p/2 ~ 1 E 



^>Y T | 2 ) p/2 



Denoting 



+ 2 3 P /2- 2]E r F ^(SYM^As*' 2 ] +2 3 P /2-2 E r f T e ^(5Y s ,5Z s dW s ) " 

L JO ^ ^ JO 



k=l 



(18) 



we apply the BDG inequality with the constant 

c-:=(^r V2 (^r>o, 

(see T heorem 3.9.1 from iKhoshnevisanl (l2002h and solution to Problem 3.29, p. 231, in lKaratzas and Shrevd 
()1995l )) and Young's inequality with some constant 72 > and obtain 



E 



eP 8 (5r„5Z.dW 8 ) 



p/2 



< E 



3=1 



e^SY' dNV 







< m p ' 2 zZ E [\J e ^ 5Y s dN l\ P/2 } < C*mP' 2 e 2 ^ s \5Yj\ 2 d{N j ) 



\p/4 



m 



( sup e^F/| 2 ) p/4 ( / e^ s d<i\P) s ) 

0<t<T JO 



p/4 



< C*m p / 2 (j 2 E 

< C*m p / 2 ( 7 2||^|| J . 



sup e^\5Y t j \ 2 ) p/2 

L 0<t<T 



+ - E 

72 



m r T 





p/2 



d(A^') s ) p/2 



< C '™ P/2+1 + ^ ll*)<^ {72||^||^ } 



/3 7 2 "•/»/ ^ 72 

where by (I14j) we have that C*m p / 2+1 = d p / 2 - With the particular choice of 

\ P/2 



(19) 



72 



:= 2 3 *>/ 2 - 1 d p/2 D- p/2 = l^/^d 



P/2 



7 



7 — dL 



>0, 



plugging (fT9j) into (fT8|) yields 



! _ ^ f /2 _ 2^ 
7 72 



<2^ 2 ~ l E (e^\5Y T \ 2 ) p/ + 2 3p / 2 " 2 E / e ^F S! 5 2 / s }ds + 2 3 f/ 2 - 2 d p/272 ||^||^, 



p/2- 
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which implies the claim. 
Step 2: We claim that 

Ds\\SYf SP <(^ p 
°p p — 2 



y/2 |( 2 P-2 + 2 3p/2-2(_^|_)f/2) ]E J(^r| (5 y T |2^/ 



_|_ ^ 2 3 P/2-2 + 2 5p/2- 



"7 — aL 



7 — aL' 
) p/2 ) E 



e f3s \(5Y s ,5 2 f s )\ds 



p/2 



holds for 



D 3 := 1 - 2 4 ?- 4 4 /2 (^) p/2 (^-) p/2 J D 2 " p / 2 - (^ r )^ ( ^)W2 2P - 2 _ 
3 p/2V p _ 2 ^ ^ 1 -6iL> 2 V 7 V p-2 y 



(20) 



(21) 



Note that the choice of K, T and a has been such that -D3 > is satisfied. To prove (|20p . we go 
back to (|16p . where we take the conditional expectation with respect to J-j, then the supremum 
over t € [0, T], raise to the power p/2 and finally apply Doob's inequality to obtain 



E 



sup (e^\6Y t \ 2 ) p/2 

0<t<T 



< E 



sup ( E 

0<t<T 



'e? r \5Y T \* + 2 [ T e^ s \(5Y s ,6 2 f s )\ds 
Jo 



< 



K p-2> 



+ / T ^(| 5 y s |2 + |^|2 )ds |^ 
Jo 7 



p/2 



y/2 ^ 2 P- 2 E[{e pT \6Y T \ 2 ) p/2 ] +2 3 ^ 2 - 2 e[( / \ (5Y S , 5 2 f s ) \ds) p/2 



+ 2 P_2 E 



— e^|«5Y s | 2 d S ) p/2 
7 



+ 2 P_2 E 



£ ^L e Ps\ 6 Z s \ 2 d S y /2 ]}. (22) 



Note that we made use of the fact that for a,b,c,d S M. and p > 2, we have 

\a + 2b + c + d\ p/2 < 2 p l 2 ~ 1 (\a + 26| p / 2 + |c + d| p/2 ) 

< 2 p ~ 2 |a| p/2 + 2 3p / 2 ~ 2 \b\ p/2 + 2 p - 2 |c| p/2 + 2 p - 2 \d\ p/2 . 

Plugging ([17)1 into (f22|) . we get 



ll^ll^ < (^-) p/2 ^~ 2 E (e^\SY T \Y /2 +2^ 2 ~ 2 E ( / e^|<«5Y s , «fe/ s )|d S ) 



y p-2' 



J T \ZV„\2\ p / 2 



\P/2 



+ 2P -2^y/^ §Y ^ + ^ty/2 D -i x 2f-2{ 2 f/ 2 E[(e^ T |«5y T | 2 ) p/2 ] 

11^1%} 



2 3p/2-l E 



< (^)P/ 2 J ( 2 P-2 + 2 3 P /2 



e^\(5Y s ,5 2 f s )\ds) 



+ 2^ 2 d 2 p/2 D^ 2 



K p-2 



-21 & L \p/2 



6lL 



1 7 — aL ' 
p/2\ 



) P/Z )E (e? T \6Y T \ 2 ) p/2 



+ (2 3p/2- 2 + 2 5 P /2-3 ( _^_ )P /2 )]E ^ y e ^|( 5 y s)52/s )| ds ) 
+ (2- 2 (^Tf/ 2 + 2^(^ \, 



\P/2 
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from which the estimate (|20p follows. 

Step 3: At this stage, estimating E[f Jq e^ s \(5Y s , 52/s)|ds 
consequence of f|20[) : Young's inequality combined with the 5?-norm yields 



p/2 

v 

0' 



will yield (|T5|) . This itself is a 



E[(^V s |(<5Y s ,5 2 / s )|d S W2 



< E 



T 



^ S \SY S \ \S 2 f s \ds) 



p/2 



<7 3 py||^ + -E 



J /3 73 



T 



e 2 



which in conjunction with the particular choice 

1 n fP - 2 \p/ 2 



(7 - aLfl 2 



p ' 2 3 P/ 2 - 2 ( 7 -aL)p/ 2 + 2 5 P/ 2 - 3 (dL)P/ 2 



> 0. 



Estimate (1201) now leads to 



^3py||^<(^) p/2 {(2^ 2 + 2*/ 2 



aL 



'p-2> 

+ ^f,p/2-2 + 2 5p/2-3 ( 



7 — dL 
5£ \p/2\ _i 



E[(e^ T |<5Y 



|2\P/ 2 1 



'7 — aL 

Notice that we trivially have \\5Y\\p < T p / 2 II 6Y\\*L P so that 

ft a Oo 



11^115* + ll^ll^ < C l p E[(e^ T \5Y T \ 2 ) p/2 ]+C 2 p E[( T e^|5 2 / s |d S ) p ] , 



where the constants Cp~ and C 2 are defined as 



Cl : =2(l+TP/ 2 ) J D 3 - 1 (--^) p/2 (2P- 2 + 2 3 f/ 2 -^ " L ^ 



C 2 := 2(1 + T p / 2 )D. 



"p-2> 
-1/ g \p/ 2 



'7 — a/7 



^ 2 3p/2-2 + 2 5p/2- 



'7 — <5L 



Moreover, it follows from (fT7j) . ([23]) and ([25]) that 

ll^l^ < C 3 E[( e ^Y T | 2 ) p/2 ] + C p 4 E[( / ef ^ 2 /.|d S ) p ] 



where the constants C 3 and Cp are defined as 



C 3 :=2D. 



-If P \P/2 D -p/2 



"p-2' 



+ (2 3 ^ 2 d 2 /2 7J- p/2 + 2 3 ^ 2 - 1 73 ) {2 p ~ 2 + 2 3p / 2 " 2 (^4 7 ) P/2 ) 



C p 4 :=2/J 3 - 1 (^-) p / 2 7J 2 ^ 2 [(2 3 ^ 2 ( i 2 + 



^ 2 3p/2-2 + 2 5p/2- 



"7 — aL 



))P/ 2 7 -l + 2 3p/2-l 



'7 — aL' 



73 



(23) 



(24) 



h^ml e^\S 2 f s \ds) p ]}. (25) 
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(recall that 73 is defined by (|24p ). From the above inequalities we obtain (I15p . where the positive 
constant C p is given by 

C p := max {C l p + C|, C 2 + C*}. (26) 

□ 

Remark 2.11. Notice that none of the constants C p , Ct and Di (i € {1, • • • ,4}) depend on the 
terminal condition or /(-,0, 0). The only problem related data they do depend on are: K, T, a and 
m. 

Remark 2.12. In the previous proof it is clear that our choices for the constants 72 and 73 do not 
lead to the most general statement of Proposition \2.9l They were chosen in this way to avoid a 
more complex statement, i.e. the constant C p given in (|26h would then depend on 72 and 73 and 
jointly with (|13p we would also have the condition D3 > 0. The conditions of Theorem \2.14\ below 
depend on the smallness of C p as given by (I26p . The particular choices for 72 and 73 lead to simpler 
expressions in our statements. 



Moment estimates - part II 



As a by-product of the two previous propositions we obtain a result on the moment estimates for 
the solution of BSDE gj). 

Corollary 2.13 (Moment estimates). Let p > 2 and (3 > 0. Let (Y, Z) £ x T-CZ be the solution 
of the delay BSDE with terminal condition £ and generator f satisfying (H0)-(H4). For K,T,a 
small enough, there exists a constant C p (which, like in Propositions \2. 7| and \2. 9i depends on several 
constants that can be suitably chosen) such that 



\Y\\ P + 

" 6 /3 



Y\ 



H ;1 



+ \\z\ 



Hp 



< cJe 



3 /3T 



\Yt\ 2 ) 



+ E 



|/M,o)| 2 d s ) p ]} 



The existence and uniqueness result 

The moment and a priori estimates in Delong and Imkeller ( 2010a! ) are tailor-made for a Picard 
iteration procedure in TL 2 x H 2 . To make such a technique work in general L p -spaces we needed to 
state a priori estimates in the form of Proposition 12.71 and Proposition 12.91 In view of those results 
one can naturally expect a compatibility condition on K, T and a more complicated than that of 
Theorem I2J2] for a solution to exist. 



With estimate (|15p at hand, we now proceed to show th e existence and uniqu e ness of solutions to 
T in S p a x TL P R for p > 2. For p = 2, Theorem 2.1 from iDelong and Imkeller] (j2010al ) (recalled in 



our Theorem I2.2p yields a sufficient condition which guarantees the standard Picard iteration to 
converge and proves the existence and uniqueness of solutions to ([!]). We will show in the following 
result that for p > 2, the convergence of the same Picard iteration is retained. What is needed to 
achieve this goal is to put up some extra effort to show that the Picard iterates (Y n , Z n ) satisfy the 
corresponding So, %^-integrability properties. 
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Theorem 2.14. Let p > 2 and assume that (H0)-(H4) hold. Let K or T or a be small enough 
such that for some /3,7 > the conditions of Proposition \2.9\ are satisfied. If further K or T or a 
are small enough such that we have 



2 p l 2 ~ 1 C v \LT J e-P s p{ds)) max{l,T p / 2 } < 1, for p € {a y , a z }, 



(27) 



where C p = C p (/3, 7, a, L, T, m) > is given by (|26p . a is given by ([7]) and L = Ka, then the BSDE 
@ admits a unique solution (Y, Z) in x H^. 

Remark 2.15. Note that, by definition of the constant C p , condition (j27j) is satisfied if either T or 
K or a is small enouqh since lim C„ < +00 which in turn implies 

lim CJaKT) p/2 = 0. 

TKa->0 

Proof of Theorem \2.14\ Let p > 2. Throughout let t € [0, T]. The proof is based on the standard 
Picard iteration: we initialize by Y° = and Z° = and define recursively 



Y t n+1 = i+ f f{s,T n (s))ds- f Z? +1 dW s , 0<t<T. 
Jt Jt 



(28) 



with F n {s) = ( f® T Y s n +v a y {dv), f° T Z™ +V a z {dv)) for s G [0, T] and n G N. In the following, let 
C > denote some generic constant which may vary from line to line but is always independent of 
n G N. We proceed by induction, where the existence of (Y ,Z ) € x T-L^ follows from classic 
stochastic analysis arguments. For n > 1, assume that (Y n ,Z n ) 6S^x solves the BSDE (]28p 
and we now prove that (|28p has a unique solution (Y n+1 , Z n+l ) £5^x H 1 ^. Note that due to 

- r T 

E 



\f(s,T n (s))\ds) 

■T 







< E 

< 2 13 - 1 E 

< 2 P ~ 1 E 
+ L p l 2 T p l 2 



T 



|/( s ,o,o)|ds+ / |/( s ,r n ( s ))-/( s ,o,o)|d s 



|/( s ,o,o)|d s J + (r/ |/( s ,r"( s ))-/( s ,o,o)| 2 d s J 



xp/2' 



|/M,0)|d S ) J 

T /-0 



y s +J 2 a y (dt;)ds + 



T r 



T 



Zn |2 



(dv)dsj 



-1P/2 



r 



T 



\Y?\ 2 ds + 



\z^d S y /2 - 



< 2 P ^E ( / |/(a, 0, 0)\ds) p + {aKT) p/2 

< 2 p - x E[( f |/(s,0,0)|ds) p j +2 p ' 2 - 1 {2aKTyl 2 {T p l 2 \\Y n \\ P sP + ||^% P ) < 00, 



(29) 



the martingale representation yields a uniquely determined process Z n+l € Hq such that 



E 



T 



£+ / /(s,r n ( S ))d S |7i =E[£+ / /(s,r»)ds] + / Z? +1 dW s , for anytG [0,T]. 
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We then define Y n+l to be a continuous version of Y t n+1 = E[£ + J. T /(s, T n (s))ds\J : 't]- Let us first 



show that Y n+1 G 5^: 



\\y 



n+li|P 



\ p - F 

cP — ^ 


' sup |y t n+1 | p 


< E 


^0 


L tG[0,T] 





sup (E[|£|+ [ T \f(s,T n (s))\ds\T t ]Y 

= [0,T1 v JO 7 



< 



te[o,T] 



,p-l 
< 2 v ~ l 



E 



/ T |/(,,r™(,))|d s 



E 



|/(,,r"( s ))|d s ) J 



< oo, 



— 1 

where the last inequality follows from the fact that £ € -L p and (|29|). This proves that y n+1 g 5q 



Since all II • llcp-norms are equivalent it follows that Y n+1 6 <S£. To see that Z n+1 G recall that 



Ito's formula applied to e^ t \Yp ,+1 \ 2 yields 



e? t \Y t n+1 \ 2 + [ (3e^ s \Y s n+1 \ 2 ds + [ e^ s \Z^ +1 \ 2 ds 
Jt Jt 

= e /3T |e| 2 + [ T 2e^{Y s n+ \f(s,T n (s)))ds- F 2e?'Q?+\ Z^dWs). 
Jt Jt 

In the above drop the two Y terms in the LHS of the equation, take t = 0, apply absolute values to 
both sides and then raise to power p/2. It follows that 



f e? s \Z2 +1 \ 2 d S ) p/2 
Jo 

<(> T |£| 2 + ^ 2e? s \Y s n+1 \ \f(s,T n (s))\ds+\J 2e? s {Y^\ Zf +1 dW s ) \ 

< 2^ 2 - 1 ( e / 3T |£| 2 ) p/2 + 2f- 2 ( F 2e? s \Y s n+1 \ \f(s,T n (s))\d S ) p/2 

Jo 

+ 2 3 P /2-2| / e ^(y s n+ \Z^ +1 dW s )\ p/2 . (30) 
Jo 



On the one hand, we have 

,-T 
H) 



E 



< E 



2e^ s \Y s n+L \ \f(s,T n (s))\ds 

T 



p/2 



T \ p/2 

2e? s \Y? +l \ |/(*, !"(*))- /(a, 0,0)|d« + / 2e? s \Y s n+1 \ |/( S ,0,0)|ds 



o 



(31) 



< c{\\Y n+ % $ e! s |/( S ,0,0)|d S ) p ] + \\Y n \\ P s , + \\Z n \\ P H$ ] < oo, 

where we have used the Lipschitz condition of / combined with calculations similar to those of ([29 
and 

[ T 2e^ s \Y s n+l \ |/( S) 0,0)|d S < sup e^|y™ +1 | 2 + ( / T e^|/( S ,0,0)|d S ) 2 . 

Jo 0<t<T v Jo 7 
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On the other hand, by the same arguments as in (]19p we find the following estimate 



E 



T 



e Ps(yn+\Z2 +1 dW s )\ 



p/2 



< d p /2 \ K \\Y 



rn+1 up 



l|| 7 n+l||P \ 
~\\ Zj W-uP (■ 



(32) 



where the last line the constant k > appear due to Young's inequality. Now choosing k > such 



that 1 - 2 2p ~ 2 dp^K' 1 > 0, it follows from (EUJ), (ED and (El that 



1 



Hg 



, ||vn+l||P 



-)ll^ n+1 ll^ <c{E[( e ^ T |£|y /2 ] 

+ E[(^ T |/M,0)|ds) p ] + ||Y»||^ + 



< OO. 



This proves that Z n+1 6 

In the next step, we prove that the sequence (Y n , Z n ) converges in x H^. Under the current 
assumptions one is able to apply a priori estimate (|15p to obtain 

||y n +i _ y n ||P _)_ \\Z n ~^~^ — Z n \\ p 



<C P E 



< C n T p/2 E 



e* s f(s,T n (s))-f(s,r n -\s)) ds 



( / e ^|/( s ,r-( s )) -/( S ,r— Hs))! 2 ^) 



\P/2 



In analogy to the calculation carried out in Equation (2.7) in Delong and Imkeller ( 2010a! ) [Proof of 
Theorem 2.1], it is easy to see that we have 



|Y 



n+l vn\\p 



<C p T p / 2 E (Lmax{ e'^ayids), J 



T 



^~-(ds)} 



e ^ a 



x (T sup e pt |Y t n - Y t n ~ l \ 2 + / - Z r s 

te[o,T] Jo 

r0 



T 



s f}s\yn yn— 1|2 



ds 



< C p T p ' 2 2 p ' 2 - 1 (Lmzx{J e^c^ds), J e~ ps a z {ds) 



p/2- 
p/2 



f T p/2prn _yn-l||P + ii^n _ gn- 



l\\P 

Hg 



/ /"0 \ p/2 

< C p 2-pI 2 - 1 (lT max {J e^ s a y (ds), J e^a^ds)}] max{l,T p / 2 } 

x ( \\Y n - Y n_1 ||^p + ||Z n - Z n - l \\ P uP) 

Hence, by ()27|) . the standard fixed point argument yields that (Y n , Z n ) converges in x which 
finishes the proof. □ 
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3 Decoupled FBSDE with time delayed generators 



The objective of this section is to extend the results from belong and Imkeller ( 2010al lbl) to the case 



of decoupled forward-backward stochastic differential equations. For measurable functions b, a, g, f, 
specified in more detail below, we study the time delayed FBSDE 



x + I b{s,X x s )ds + [ a(s,Xf)dW a , x£R d , (33) 
Jo Jo 

Y t x = g(X%)+[ f(s,Q x (s))ds- [ Z*dW s , 0<t<T, (34) 
Jt Jt 



where for t € [0, T], we write 

Q x (t) = ({X* ■ a x )(t), {Y* ■ a y )(t), (Z x ■ a z )(t)) 

= [J T Xf +v a x (dv), f^Y t \ v a y (dv), J^Zf +v a z (dv))^ (35) 

with given deterministic finite measures a x ,a y and a z supported on [— T, 0). The coefficients 
b,a,g,f appearing in (j33p ~()34 |) are assumed to satisfy certain smoothness and integrability condi- 
tions such that the backward equation (|34p falls back into the setting of (H0)-(H4) from Section 
12.11 More precisely, we assume the following to hold: 

(F0) a x , a y ,a z are three non-random, finitely valued measures supported on [— T, 0); 

(Fl) g : R d — > R m is continuous differentiable with uniformly bounded first order derivatives, i.e. 
there exists K' > such that |Vg| < K'; 

(F2) / : [0,T] xM^x R m x R mxd R m is continuously differentiable with uniforml y bou nded 
derivatives, i.e. there exists a constant K > such tha10 \V x f\, \V y f\, |V Z /| < yJ~K/?> holds 
uniformly in all variables; / satisfies a uniform Lipschitz condition with Lipschitz constant 

(F3) b : [0, T) x R d -)• R d and a : [0, T] x R d -> R dxd are continuously differentiable functions with 
bounded derivatives; |6(-,0)| and |<r(-,0)| are uniformly bounded; a is elliptic; 

(F4) (/ T |/(s,0,0,0)| 2 d S ) p/2 < oo for p > 2; 

(F5) /(i,v,-)l(-oc,o)(i) = 0; 

Condition (F3) is a standard assumption which guarantees the existence and uniqueness of the 
solution of SDE (|33p . Furthermore, condition (F2) implies that the generator is uniformly Lipschitz 
continuous in (a;, y, z) G R d x R m x M mxd . In analogy to conditions (H2) and (H2') from section [2~Tl 
let us write down the following implication of the Lipschitz condition (F2): with the constant K > 



1 We remark that this bound is taken over the corresponding Euclidean norm of the derivatives matrix/tensor. 
To avoid possible confusion when using tensors one can always interpret / in the variable z € R mxd as taking 
not a matrix but a sequence of d-dimensional vectors Zi € R d (i £ {1, ■ ■ ■ ,m}). The condition would then read 
Efci Nnf \ < V^/S where / : [0, T] x R d x R m x R d > 
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chosen above, for any t £ [0, T] and any sufficiently integrable vector or matrix valued processes 
u, u', y, y' and z, z' it holds that 



(F2') |/(t, (u • a x )(t), (y ■ a y )(t), (z ■ a z ){t)) - f(t, {v! ■ a x )(t), (y' ■ a y )(t), (z' ■ a z ){t))^ 
<K(\(u-a x )(t)-(u'-a x )(t)\ 2 

+ \{V <x y )(t) - {y' ■ a y )(t)\ 2 + \(z- a z )(t) - (z 1 ■ a z )(t)\ 2 ) 
< Ka x ([-T,0]){(x - x'f ■ a x )(t) + l((( V - y'f ■ a y )(t) + {(z - z'f ■ a z )(t)) 

where L := Ka with a defined in ([3]). For a fixed x £ M. d , the existence and uniqueness of solutions 
to the backward equation (jM|) in <Sj x 'Ho is guaranteed under the assumptions (F0)-(F5) together 
with the compatibility criterion from Theorem 12.21 on the terminal time and the Lipschitz constant 
L = Ka, i.e. 

1 f'° 

(8T + -)L / e" /3s p(ds)max{l,T} < 1, for p G {a y ,a z }. 

P J-T 

To extend the result to 5^ x T-fc for p > 2, one only needs to replace the condition above by the 
compatibility condition from Theorem 12.141 

/ /"O \ p/2 

T > I 2 ~ 1 C P [LT J e- /3s p(ds)J max{l,T p / 2 } < 1, for p G {a y ,a z }. 

Throughout this section, given p > 2, we will assume that for every x G M d , the FBSDE (|33|) -(|34|) 
admits a unique solution (X x , Y x , Z x ) G S^W 1 ) x S p ^(R m ) x V p f} {W nxd ) for all q > 2. 

3.1 Gateaux and Norm differentiability 

In this section we investigate the variational differentiability of the solution (X x , Y x , Z x ) of the time 
delayed FBSDE (}3"3")) - (f3"4"|) with respect to the Euclidean parameter x G M. d , i. e. with r espec t to the 
initial condition of the forward diffusion. By a well known result (see e.g. iProtterl ()2005l )). (F3) 
implies that the forward component X x is differentiable with respect to the parameter x G It 
is natural to pose the question whether this smoothness is carried over to (Y x , Z x ) in the setting of 
FBSDE with time delayed generators. In all this section we fix h an element of M. d \ {0}. Our goal 
is to show that the variational equations of (|33j) -(|34|) are given by 



VXfh = h+ f Vb(s,X x )VX x hds+ I Va(s,X x )VX s hdW s , (36) 
Jo Jo 

VY t x h = Vg(X%)VX%h - [ VZ*hdW s + [ ((V f)(s,G x (s)) , (VQ x h)(s))ds, (37) 

Jt Jt 

where the notation VX 1 (respectively VY X and V ' Z x ) denote the Gateaux derivatives of X x (re- 
spectively Y x and Z x ) in the direction h and (VQ x h)(t) is to be understood in the same fashion as 
in (f35|) . i.e. 

(V& x h){t) = ({VX x h-u x ){t),{VY x h-a y ){t),{VZ x h-u z ){t)), t G [0, T] . (38) 
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Note that (F3) implies that (i36|) admits a unique solution in Sp for every (3 > and p > 2. Let 
(X,y,Z) and VI/i solve (j53>(jMl) and ([36]) respectively and let 9^ be as defined by flM}. Now 
consider the BSDE with the linear time delayed generator for t £ [0, T] 

P t h = Vg{X?t)VX?th- f Q s hdW s + [ P(s,(Ph ■ a y )(s),(Qh ■ a z )(s))ds, (39) 

Jt Jt 

where F : U x [0, T] x R m x M mxd R m , F(t,p,q) = ((V/) (t, 9 x (t)) , ((VI^-a*)(t),p,?)). 
The next corollary states, using Theorem 1 2 . 2 1 and Proposition 12.71 a result concerning the existence 
and uniqueness of solution to (|39p . This solution process will then serve as the natural candidate 
(in some sense) for V ' r Y x h and V T Z x h, solution to ([37. 



Corollary 3.1. Let p > 2, h £ R d \ {0} and /3 > 0. A ssume that (F0)-(F5) are satisfied and let 
L > be as in (F2'). If p > 2 assume that T, K, a are chosen like in Proposition \2.9\ and satisfy 
in addition 



2 p / 2 ~ l C p [LT J e-P s p(ds)j max{l,T p / 2 } < 1, for p € {a y ,a z }, 



(40) 



If p = 2 assume T, K , a are chosen such that the conditions of Theorem \2.2\ and of Proposition 2.1 
hold. Then for every fixed x in M. d , BSDE (|34p has a unique solution (Y, Z) € x 7-L P p and BSDE 
(I39p has a unique solution (Ph,Qh) £ S 1 ^ x ~H^. 

Proof. Given the known properties of X and VJ (and hence of Vl/i) it is easy to see that £ = 
Vg{X*)VX*h and F(-,0,0) satisfy conditions (HI), (H3) and (H4). We recall Remark EH] to say 
that the several compatibility conditions (|40p as well as the conditions in Proposition 12.91 depend 
only on the Lipschitz constant K of (F2), the delay measures a y , a z , T and the dimension of the 
equations. 

From the definition of F and using the bounds of the (spatial) derivatives of / assumed in (F2) it 
is clear that F satisfies a standard Lipschitz condition (in the spatial variables). In particular, take 
p,p' € M m and| q,q' G R mxd , then via Minkowski's and Cauchy-Schwarz inequalities along with 
(F2) we have 

\F(t,p,q)-F(t,p',q')\ < \((V y f){t, &°(t)) , (p - p')}\ + |((V,/)(t, &°(t)), (q - q'))\ 

< \(V y f)\\p-p'\ + \(V z f)\\q-q'\ < y^K/3( \p - p'\ + \q - q'\ ). 

And hence F satisfies exactly the same Lipschitz condition as /. Furthermore, the delay measures 
appearing in F are exactly the same ones as those that appear in /. We can thus conclude that the 
Lipschitz constant, the delay measures, terminal time T and dimensions for / and F are the same. 
Under this corollary's assumptions, the conditions of Theorem 12.141 are satisfied for both BSDE 
(|34l) and (d|. The existence of a unique solution (Y, Z) and (Ph, Qh) in x Hp of ([34]) and (J39l) 
(respectively) follows from Theorem 12.141 (and Theorem I2.2p . □ 

The solution of BSDE (I39p serves now as the natural candidate for the variational derivatives of 
(Y, Z) solution of (|37p . If one shows that (VY x h, V ' Z x h) exist in some sense then by the uniqueness 



2 Or a sequence of q[ £ R m with i 6 {1, • • • , d} as we saw in page 1201 s footnote. 
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of the solution of ([39]), the solutions to ([37]) and ([39]) must coincide, i.e. (VY x h, VZ x h) = (Ph, Qh) 
holds almost surely. 

For the rest of the section, we assume that all assumptions ensuring the existence and uniqueness 
of the variational equations ()36[) - (|3T[) are fulfilled, i.e. we assume that the assumptions of Corollary 
13.11 hold. In our next result we show the mapping x i-> (Y X ,Z X ) is differentiable in an adequate 
sense. 



Proposition 3.2. Take p > 2 and assume the conditions of Corollary (|3.ip hold. Then for any 
x G R d the solution (X X ,Y X ,Z X ) of the FBSDE is norm- differentiable in the following 

sense: 



lim 



Y 



x+eh 



Y 1 



VY x h 



CP 



lim 



7x+eh 



Z 3 



VZ x h 



0, Vh£R d \ {0}, 



nil' 
H ;:l 



where (S/Y x h, VZ x h) is the unique solution of the BSDE 



A 



s,X 



VY x h = Vg{X^)VX^h- [ VZ x hdW s + [ <(V/) (s, Q x (s)) , (V@ x h)(s))ds, 

Jt Jt 

with Q x and VQ X defined by (|35p and (|38p respectively. 

Proof. Let x G M, d , t G [0,T] and e > 0. We use the following notations 

J V x f(s, (X x ■ a x )(s) + 6((X x+£h - X x ) ■ a x )(s), 

(Y x ^ h -a y )(s),(Z x+£h -a z )(s))d9, 
j\ y f(s,(X x -a x )(s), 

(Y x ■ a y )(s) + 9({Y M - Y x ) ■ a y ) (s), (Z*+* h ■ a z )(s))d6, 
J V z f[s,(X x -a x )(s), 

(Y x ■ a y )(s), (Z x ■ a z ){s) + 8{{Z x+eh - Z x ) ■ a z )(s))d6. 



y 



(41) 



A 



s,Z 



We remark that although the processes A depends on e and x, for the sake of notational simplicity 
we do not write this dependence explicitly. We remark also that by assumption (F2) the processes 
1^.-,* I < \J K/2> for * = X,y, Z, in particular they are uniformly bounded in x and e. 
We denote by (Ph, Qh the solution of the BSDE ([39]) which coincides with (VYh, VZh). We define 
the auxiliary processes f := (g(X x+sh ) - g(X x ))/e - Vg(X x )VX x h, 

i^x+eh ~y x yx+eh yx vx+eh vx 

U := Ph, V := Qh, and X := VX x h. (42) 

e e e 

Notice that from Assumption (F2) and the standard SDE theory we have that X is well defined 
and for any b > and p > 2. We now claim and prove that 



lim ||?7||pp = lim ||V||?.p = 0, for arbitrary x G 

rt 



£->0 



5 /3 £->0 
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This result obviously proves the norm differentiability. To start with, we have 

T f(s,e^ h (s))-f(s,@-(s)) 



u t = t- 



-da 



£ 



r 



;(V/)(s, O x (s)), ((VX x h ■ a x )(s), (Ph ■ a y )(s), (Qh ■ a z )(s)))ds 



VAW X . 



By construction the above equation is well defined, since for any x and e all the involved processes are 
known a priori to exist and have the convenient integrability properties. The format of the above 
dynamics is still not convenient for our computations so we transform it into the more familiar 
dynamics of a delay BSDE. Using the identity 4>{x) — 4>(y) = (x — y) J V(f)(y + 6(x — y))dO for a 
continuously differentiable function <fi : IR a — > W b (a and b being arbitrary non-zero integers), the 
previous equation leads to 



U t = £ + - / [A SiX ((X*+ £h - X x ) ■ a x )(t 
£ Jt 



+ A s , y ((Y x+£h - Y x ) ■ a y )(s) + A S)Z {{Z* +eh - Z*) ■ a z )(s)]ds 



[(Vf){s,e x (s)),{(VX x h • a x )(s), (Ph ■ a y )(s), (Qh ■ a z )(s)))ds 



T 



V,dW, 



T 



^+ / $ (s,(X ■ a x )(s),(U ■ a y )(s),(V ■ a z )(s) )ds 



T 



YAW, 



(43) 



with X given in (gSJ), $(t, x, y, z) := R t + xA t)X + yA t ,y + zA tyZ and 

R t := -((V/)(t,e^)), ((VX x h ■ a x )(t), (Ph ■ a y )(t), (Qh ■ a z )(t))) 

+ A t:X (VX x ■ a x )(t) + A t: y(Ph ■ a y )(t) + A t>z (Qh ■ a z )(t). 

We now aim at using the results of Section 2 on the family (index by e) of auxiliary delay BSDEs 
(|43p . In view of the uniform boundedness of the processes A and the linearity of the driver we 
can repeat the arguments used in the proof of Corollary 13. II to conclude that under the assumptions 
of this proposition the data of BSDE (|43p (Lipschitz constant, delay measure and terminal time) 
satisfies uniformly in e the assumptions of Corollary 13.11 as well. 

Applying the a priori estimate of Proposition 12.91 or the moment estimate from Corollary 12.131 to 
the BSDE (|43p and taking into account that $ satisfies (F2), we get 



\U\\ p s v + \\V\\ P H v <c p {\ 



<C\E 



(e? T \t\ 2 y' 2 

(e0T|£|2)P/2 



+ E 



+ \\X\\l P + E 



>(s,(X-a x )(s),0,0)\d S ) p }} 
j\P s \R s \ds) p }}, (44) 



for some constant C > (where we have used that A. t x is uniformly bounded). We proceed to 
compute the limit of each term on the right hand side of (|44p as e goes to zero. 

We first deal with the second term of the right hand side of (|44p . Define 



0"t := 



Va (t, X x + 9(X x+£h - X?))d6 and b t :-- 



Vb(t, X x + 6(X x+£h - X x ))d9. 
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Note that X € S p for any p > 2 (see (|3Sjl ) and solves the linear SDE 



X t = J 4 + / [ffj,]d^+ / [b s X s ]ds, 



(45) 



where J is given by 



Jt :-- 



o 



~yX*h(a s - Va(s,X*))]dW s + / [ VX*h(b s - Vb(s,X*)) ]ds. 



o 



Given the known properties of VX and the fact that b, a, V6, and Vcr are uniformly bounded we 



have that J G <Sa for any p >2. Indeed, Doob's inequality leads to 



E 



sup I / [VX x s h(a s - Va{s,X*))]dW s \ 2 Y /2 ] < C\\ VX x h(a - Vtr(-, X x )) II* < oo. 
te[o,T] Jo 7 J 



Moreover, note that by Lebesgue's dominated convergence theorem 

lim \\VX x h{a - V<t(-, ||^ p = 0. 

Similarly, using Jensen's inequality, the finite variation part of J is an element of <Sq(R) and 

lim II Jll qp = 0. 

Now we derive the following estimate for X in terms of the norm of J 



\X\\ s p <CE[ sup |X t | p ] < C \\J\\ s p 
13 te[o,T] 



(46) 



which will show that lim e ^o ll^ll<s p = 0- Indeed equation (|45p implies that: 



E[ sup \X r \ p ] <CE 

0<r<t 



sup 

0<r<t 



0<r<t JO 



J r \ p + sup I / [a s X s ]dW s \ p + sup I / [b s X s ]ds\ p 



0<r<t JO 



Applying Burkholder-Davis-Gundy inequality to the second term in the right hand side, we get: 



0<r<t 



-0<r<t 



E[sup \X r \ p ]<CE sup |J r | p + | / |a s X s | 2 ds| p/2 + sup I / [b s X s ]ds\ 



0<r<t Jo 



Jensen's inequality and the fact that a and b are bounded imply that: 



E[ sup \X r \ p ] <CE sup |J r | p + / \X s \ p ds 

0<r<t L 0<r<i 



hence 



E[ sup \X r \ p ] <c{e[svcp \J r \ p ]+ [ E[ sup |X r | p ]ds). 
0<r<t ^ 0<r<t Jo 0<r<s > 



Gronwall's lemma finally entails estimate (|46p and thus lim e ^.o H-^lls? = 0. 
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Let us consider the terminal condition term in (I44|) . Denoting 



g:= ! Vg{X% + e(X*+* h - X£))d9, 
Jo 



it holds that 



E 



(e /3T|£|2 )P /2 



yX+eh _ yx 

e /3Tp/2 1 1 - / T T 



e 

yx+eh yi 

< C JH A T T 



VX|/i) + (g- Vg{X^))VX^h 



xu\\P 
LP 



'{II 



VX^h ' + \\\VX£h\ \g-Vg(X%)\ 



X \ I MP 

I LP 



< C{||X||^ + || |VXf/*| \g - Vg(X%)\ \\ P LP }} ^ 0, 

where we have used Lebesgue's dominated convergence theorem for the second summand and the 
estimate obtained above on the norm of X for the first one. 



Now, let us consider the last term on the right hand side of (|44p . We have that 



E 



[(/Vi R ,i ds ) p ; 






< CE 











+ CE 
+ CE 



e^ s | (A SjX ~ V x f{s, @ x (s))) (VX x h ■ a x )(s)\ds 
( £ e?' | {A.y - V y f(s, @ x (s))) (Ph ■ a y )(s)\ ds 



A 



s,Z 



V z f{s,@ x (s))) {Qh-a z ){s)\dt 



Standard arguments yield (note that e > is implicitly contained in A^Xi see (|4ip ) 
A t ,x — > V x /(t, Q x (t)) as £ -> in probability, for dt-a.a. t € [0, T]. 
Moreover, Proposition 12.91 and the previous calculations show that 



\y 



x+eh vx\\P 



Y X K P + \\Z X 



x+eh 7%\\V 



\ijp 



< C {ef T P\\g{X x+£h ) - g(X x )\\ p LP + \\X x+eh - X x f P } — ► 



fi E-tO 



for some positive constant C. This implies for dt-a.a. t € [0, T] 



x+eh 



-+ Y x , Zf +eh -+ Zf, as £ -+ in probability. 



■t ^ ± t i 

Since Vy/, V 2 / are continuous, it follows that for di-a.a. t £ [0, T] 

A t y — ► V y f(t, ® x {t)), as £ — > in probability, 
A t .z — > V z f(t, @ x {t)) , as £ — > in probability. 

Thus, using Lemma 11.11 and the fact that P and Q are square integrable, Lebesgue's dominated 
convergence theorem (which also holds, if almost sure convergence is replaced by convergence in 
probability, cf. IShirvaevI (|1995l ). remark on p. 258) yields lim^o E [( /c f e^ s \R s \ds) p ] = 0. Now 
yields that 

lim{||tf||* P + ||y||L,} =0, 



which proves the claim. 



□ 
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3.2 Strong differentiability 

All previous assumptions on existence and uniqueness remain in force. In this section, we concen- 
trate on the smoothness properties of the paths associated to the processes (Y X ,Z X ). We assume 
throughout this section that m = 1, i.e. the delay BSDE are now one-dimensional. A first result is 
obtained in the following 

Proposition 3.3. Set m = 1 and under the assumptions of Corollary \3A\ we have for x,x' 6 R rf 

E[ sup \X X - Xf'\ q ] < C\x - x'\ q , for any q>2, 
0<t<T 



and for any p > 2 
E 



sup (c*|l?-lf| a ) 

0<i<T 



2\P/1 



+ E 



e? 8 \ZZ-Z?\ 2 d 8 ) p/i 



< C\x-x'\ p . 



Thus for every x £ 



• the mapping x t— > Y x from M. d to the space of cddldg functions equipped with the topology given 
by the uniform convergence on compacts sets is continuous V-almost surely, 

• the mapping x i— > Z x is continuous from M. d to L 2 ([0, T]) F-almost surely. 
In particular, for every x € M. d , 

• the mapping x i-> Yf from R rf to M. is continuous for all t E [0, T], ¥ '-almost surely, 

• the mapping x i-> Zf(uS) is continuous for every x £ M d and dt <S> dF-almost all (t,u). 



Proof. The estimate on the forward process is classical (see e.g. (jProtterl . 120051 . Theorem V.37 
Equation (***) p. 309)). In this proof, C > denotes a generic constant which may differ from line 
to line. We apply the a priori estimate from Proposition 12.91 and get 



E 



sup (e^\Y t x -Y t x '\ 2 ) p/2 



0<t<T 





- r T 


+ E 


( / e? s \Z x 




- Jo 



rx'\2 An \P/ 2 



{e^\g{X x T )-g{X x ^y 



2\p/2 



+ E 

< c{e 

+ E 



"( ^ el s |/( S , (X x ■ a x )(s), COO) - f{s, {X x> ■ a x )(s), ((s)) \ds) p ] } 

y T \g(X^-g(X^\Y /2 ' 
\ £ eP'\f(8, (X x ■ a x )(s)X(s)) - f(s, (X x ' • a x )(s), C(s))\ 2 ds) p/2 ] }, 

with C(-) := {(Y x ' ■ a y )(-),(Z x ' ■ a z ){-)). Using the mean value theorem and the boundedness of 
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V/ and Vg (i.e. the Lipschitz property of / and g), we deduce 



E 



sup (e^\Y t x -Y t x '\ 2 ) 

0<t<T 



2\P/2 



(e^lXZ-X 



+ E 

-x'|2-jP/2 



+ E 
+ E 



e? s \Z*-Zf \ 2 ds) p/2 
j\^\((X x -X x ')-a x )(s)\ 2 dsY /2 ]} 
ef s \X x -X x '\ 2 ds) p/2 ]} 



< C\x-x'\p, 



where the last two lines follow by applying the change of integration from (jlOp and the first claim 
of the proposition. The continuity properties of the mappings x 1-4 Y x and x i -> Z x are now 
obtained by an application of Kolmogorov's continuity criterion (see for example (jProtterl . 120051 . 
IV.7 Corollary 1)). □ 

If the generator exhibits additional regularity, it even turns out that the paths of x 1— > Y x are 
continuously differentiate. 

Theorem 3.4. Let /3 > and assume the conditions of Proposition \3.S\ can be verified for some 
p > 4. Assume moreover that all (spatial) second order partial derivatives ofb,a,g and f exist, 
are continuous and uniformly bounded. Then, for any (x,e),{x',e') G R d x (0,oo), h G R d and 
p G (2,p/2] it holds that 



E 



sup I e H 

0<t<T 



x+eh \-x 



-w-x'+e'h \sx' 



2\p/2 



j < C {\x-x'\ 2 + \e-e'\ 2 ) 



|2\P/2 



Thus V X Y X belongs to 7i p ^ and the mapping x 1-4 Y x (uj) is continuously differ entiable for all t € 
[0,T], F-almost surely. 

It is known that the existence of the partial derivatives (or even all of the directional derivatives) 
of a function does not guarantee that the function is differentiable at a point. But it is if all the 
partial derivatives of the function exist and are continuous in a neighborhood of the point, then the 
function must be differentiable at that point and is in fact of class C 1 . 

Under the assumption that m = 1 and the subsequent corollary of the Theorem in the previous 
section, we know that the all (spatial) partial derivatives of Y x exist. The main result of Theorem 
3] is the continuity of those partial derivatives. 



Proof. As in the previous proof, C > denotes a generic constant which can differ from line to 
line. Let p > 2, t G [0, T] and h G R d \ {0}. For (x,e) G M. d x (0, 00) let U x > e := y*+"*-y* ^ 



yx,e z xJreh —Z x ^x,e 



g(X*+ sh )-g(X*) 



— and X x,£ := xx+£h xx _ Using the notation from the proof 



of Proposition the pair (U x ' £ , V x ' £ ) satisfies the BSDE 

JJ x t > £ = t x > £ + ^ ( x > £ ( s ))ds - ^ V t 



dW, 



with C' E (t) := ((U x > £ ■ a y )(t), (V x ' £ ■ a z ){t)) and $(t,y,z) := (X x > £ ■ a x ){t)A x f x + yA%$ + zA x %. 
Note that the terms A x f with * = X,y,Z are given by ([41 p . 
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For whatever choice of (x,e) we emphasize that the arguments used in the proof of Corollary 13.11 
and Proposition 13.21 hold true for the above auxiliary BSDE in what the applicability of the a priori 
estimate of Proposition 12.91 is concerned. 

Let another pair (x',e') G R d x (0,oo) be given. Applying Proposition 12.91 yields 

\\U X > £ - U x '> e '\\ p s v < C p \E[(e^ T \^ £ -e'' £ '\ 2 Y /2 } e2 s \5 2 <P(s)\ds) p } }, 



with 



5 2 $(t) := (X*' e • a x )(t)J%% - (X x ' £ ■ a x ){t)A x '> 



t.X 



+ (U x > £ • a y )(t)(AQ - A X J ) + (V* > £ ■ a z )(t){A x f z - A X J ). 
Using the hypotheses on / (i.e. all partial derivatives up to order two are bounded), we find 

|da*(t)| < c{\((X x > £ - X x '' £ ') ■ a x )(t)\\A x ; x \ + \(X< £ ' ■ a x )(t)\\A x ; £ x - A?/\ 

v',e'- 

As a consequence 



+ \(V^'.a z )(t)\\A^ z -Af/\}. 



px,e _ jjx ,e ||P 



< c{\\^-Z*'>*'\\I p+ e[( / eH(X^-X x '^-a x )(s)\\A x s f x \ds" p 



+ E 
+ E 
+ E 



( / e^'\(^ - a x ){s)\\A^ - A'J \ds) 1 



T 



Hu x '' £ '-a y )(s)\\A^-A^'\dsy 



T 



e* s \(V*> £ -a^WA 



x.e a x 



A 



;,Z l ds ) J } 



<c\\\e' E -i 



x ,e \\P 
1 1 LP 



x.e yx ',£' lip II AX,£ up 



+ \\X X > £ -X 



,, A' 



+ Wx 3 



H 



I ax,e _ a: 



I \\TJX 1 ,e' ||P || AX,£ aX ,e up J_ \\\ix' ,£' \\p \\ ax,£ ax 1 ,e up 

where for each term we used the Cauchy-Schwarz inequality twice, that A 1 < e^ 1 and (jlOp . Since 
(U* > £ ,V X ' £ ) is a solution in S p x % p of a BSDE, it follows from Corollary 12.131 that the quantities 

E[( J T e^ s \uf' £ '\ 2 ds) p ] and E[( J q T e /3s |V r /' e '| 2 ds) p ] are finite and uniformly bounded in e' . By the 
assumptions on b and a, we have 



E 



V|A^'| 2 d S ) pl1 



< oo. 
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In addition, by the boundedness of V/ we have that \A X *\ and \A X * £ | are uniformly bounded (in 
their several parameters) with * = X,y, Z. Thus the estimate reduces to 



II _ 4a; ,e MP 1 || A x < s _ A 



x' ,e' up 1 



(47) 



Using the mean value theorem and the fact that the second order partial derivatives are bounded 
it holds that 



A x,e _ ax ,e I , I ax,e A \ -\- I A 



A+z 



<c{(\ 



x x+eh _ x x>+e'h\ . ax j(Q + ^yx+eh _ yx'+e>h\ . Q ^ 



+ (l^ 



x+eh yx'+e'h\ 



a z ){t) + (\X x -X x '\-a x ){t) 



+ fly* - y*'| • a y )(i) + (\Z X - Z*'\ ■ a z )(t)}. 
Plugging the right hand side of this inequality in ()47|) and using Lemma ll. II one gets 

\\u x ' £ - u x '> e 'f qP < c\u x ' £ - i x '' £ '\\ p L v + \\x x > £ - X x '' £ '\\ p 2p + \\X X - x x ' ' 



U 2 P 



+ II* 
+ \\z 



Y 1 



x+eh jg-x'+s'h ||P _|_ ||ya:+f 

x+eh ryx'+e'hup 



H 



H 



_l_ \\V X V x ' \\P J-\\7 X 7 X ' \\P I 
2p + ll r -I \\aj2p + \\^ \\ v 2p(- 
tt a rta J 



Since b, a and g are twice continuously differentiable with bounded derivatives we have the following 
estimate 

E[\C X ' £ - r V r] < C(\x - x'\ 2 + \e- e'\ 2 ) p/2 , 

which is proved for example in (jAnkirchner et al.l . 120071 . Lemma 7.4). This result combined with 
Proposition 13.31 leads to 



E 



/3t\ TT x,e _ tjx' ,e' \2\p/^ 



sup (e^U^-U, 



0<t<T 



<C(\x-x'\ 2 + \e-e' 



2NP/2 



The last claim of the theorem follows using Kolmogorov's continuity criterion (see for example 
dProtterl, \2*M, IV. 7 Corollary 1)). □ 



4 Representation formulas and path regularity 

One of the fundamental results in the setting of FBSDE concerns the relationship between the 
Malliavin and the variational (classical) derivatives of the solution process: the Malliavin derivative 
of the solution of the BSDE can be expressed as a product of the BSDE's solution variational 
derivatives (with respect to the initial parameter of the SDE) and the variational derivatives of 
the forward diffusion. This relat ionship is known to ho ld both in the standard Lipschitz generator 
setting (see Proposition 5.9 of El Karoui et al. ( 19971 )) as well as the quadratic generator case 
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(see e.g. Theorem 2.9 of llmkeller and Dos Reisl (l2010ah ) for classical BSDE without time delayed 
generators. 

In this section we show that this relationship still holds for decoupled FBSDE with time delayed 
generators. Such a result is somewhat surprising since it is normally dependent on a Markovian 
structure for the solution of the BSDE that exists for non-time delayed BSDE and which fails to 
materialize for time delayed BSDE. Imperative for this relationship to hold is the fact that the 
forward process X is Markovian along with a good behavior of the terminal condition. 

As in the previous section, whenever we consider the delay FBSDE (|33p - (f34"l) . we assume that all 
conditions to ensure the existence of a unique solution {X, Y, Z) are in force. Moreover, since for 
ft > 0, all /5-norms are equivalent, in the following we content ourselves with giving results for ft = 0. 
Recall that we assume m = 1, i.e. the delay BSDE is not vector-valued. 

Malliavin's differentiability of FBSDE with time delayed generators 



We recall Theo rem 4.1 oflDelong and Imkellerl (l2010bh . modified to our the FBSDE setting. The 



orem 



4.1 from Delong and Imkellei (|2010bl ) shows that the solutions of time delayed BSDE 



are 

Malliavin differ entiable, and as a consequence, it can be deduced that the solution of the time de- 
layed FBSDE (|33p - (|34p is also Malliavin differ entiable. Under the condition (F3) on the coefficients 
of the forward equation (|33j) . the Malliav in differentiabi lity of the forward process X is a standard 
result, see for instance Theorem 2.2.1 in iNualartl (jl995h . We denote the solution to the equations 



([33])- ([34]) by (X, Y, Z). The next result states the Malliavin differentiability of (X, Y, Z). Using the 
notation introduced in Section 3, we define for < u < t < T 



(D u e)(t) = ((D U X ■ a x )(t), (D U Y ■ a y )(t), (D U Z ■ a z ){t)) 



\ D u X t + v a x (dv), / D u Y t+v a y (dv), / D u Z t+v a z (dv ) ) . 

J-T J-T J-T J 



(48) 



We define in the canonical wa>H the space Li 2 as the space of progressively measurable pro- 
cesses, X € H 2 , that are Malliavin differentiable and normed by ||-X"||l 12 = EfJ^ |X s | 2 ds + 
/ T J T \D u X s \ 2 dsduYl\ 

Theorem 4.1. Take p = 2, m = 1 and assume the conditions of Corollary \3.1\ hold. Then (X, Y, Z) 
are Malliavin differentiable and their derivatives (DX, DY, DZ) solve uniquely in Li 2 x Li 2 x L1.2 
the following time delayed FBSDE: 

D u X t = a(u,X u )+ I V x b{s,X s )D u X s ds+ i V x a{s, X s )D u X s dW s , (49) 

J u J u 



D u Y t = Vg(X T )D u X T - C D u Z s dW s + C <(V/) (a, 6(s)) , (D u Q)(s))di 

Jt Jt 



(50) 



forO < u < t < T (zero otherwise) with and DQ given by (]35p and (j48|) respectively. Furthermore, 
{D t Y t : t G [0,T]} is a version of {Z t : t € [0,7]}. 



s See Section 2.2 of llmkeller and Dos Reisl (|2010aT) . Section 5.2 of lEl Karoui et all (|l997l ) or simply INualartl (|l995l ) 
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Proof. The results concerning the forward component are well known, s eelNualartl (119951) orllmkeller and Dos Re 
^OlOal ). The conditions of Corollary O ensure that Theorem 4.1 from lDelong and Imkellerl (|2010bl ) 
can be applied. Hence Y and Z are Malliavin differentiable. The representation of Z by the trace 
of of the Malliavin derivative of Y follows as well from the cited result. □ 

The representation formulas 

We now present the representation formulas for (I49p and (|50p which are effectively expressed in 
terms of the variational VX, VY and VZ. 

Theorem 4.2. Let the conditions of Theorem 4.1 hold. Let (X, Y, Z), (VX, VY, V Z) and {DX, DY, DZ) 
denote the solutions of FBSDE f)33[) -(|34 |) . (|36p -(j37 p and (|49p -(l50 |) respectively. Then the following 
representation formulas hold: 



D u X t = VX t (VX u )- l a(u,X u )l {u < t} , 
D u Y t = VY t {VX u y l o-{u,X u )l [u < t}) 

Z t = VY t (VX t )- 1 a(t,X t ), 
D u Z t = VZ t {VX u )' l a{t,X u )l [u < t} , 



t,u€ [0,T], dP- a.s. (51) 
t,u e [0,T], dP- a.s. 
t G [0,T], dPOdt- a.s. (52) 
t,uG [0,T], dPOdt - a.s. 



Proof. As in Theorem 14.11 we remark briefly t hat the proper ties of the forward component ar e 
well known and hence equality (I5ip holds, see iNualartl (j 19951 ) or llmkeller and Dos Reisl (|2010al ). 
Theorem 14.11 ensures that (DX, DY, DZ) is the unique solution of the time delayed FBSDE (|49p ~ 
([50]) . Throughout let t € [0, T] and u G [0,t]. We define the processes 



U Uyt = VY t {VX u )- l a{X u )l {u < t] and V u , t = VZ t {VX u )~ 1 a{X u )t {u < t} , 
and for s € [0,T], we set D u X(s) = f® T D u X s+v a x (dv), 

U u (s)= / U UtS+v a y (dv) = / S/Y S+V (S/X U ) a(u, X u )t{ u < s+v} a y (dv ), 

J-T J-T 

K(s) = J V u , s +va z (dv) = J VZs+^VX.^o-^X^t^s+^azidv), 

compare also with the notation in (pQ). Multiplying the BSDE ([37p with (VX u )~ 1 o~(u, X u ) and then 
using (|5ip we obtain for any < u < t < T dP-a.s. that 

r-T 

*dW s 



U u>t = Vg(X T )D u X T - J V u , 



+ 



{(Vf){s,e(s)),(D u X(s),U u (s),V u (s)))ds, 



where Q is given by G(-) = [(X ■ a x )(-), (Y ■ a y )(-), (Z ■ a z )(-)) (compare with (|35|) from section 
[3|). Now, Theorem 14. II states that the solution of BSDE ([3D]) is unique, hence (U, V) must coincide 
with (DY, DZ). Another way to see this would be to use the a priori estimates of Proposition 12.9 
with ([3D]) and the above BSDE. 

Formula (|52p follows easily from a combination of the representation formula for D u Y t combined 
with D t Y t = Z t , dP <g> dt-a.s. (see Theorem EI}. n 
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Implications of the representation formula 

The representation formulas in the previous theorem allow for a deeper analysis of the control 
process Z concerning its path properties. 

Theorem 4.3. Let p > 2, assume that |/(-,0, 0, 0) | is uniformly bounded and that the conditions of 
Corollary 13.11 hold. Then for p > 2, the mapping t i— >■ Zt is continuous dP-a.s. If moreover we have 
p > 2, then we also have 

Hulls' <°° forqe [2,p). 

In particular, for p > 2 we have for every s, t G [0,T] that E[\Y t - Y s \p] < C\t - s\ p / 2 and that Y 
has continuous paths. 

Proof. It is fairly easy to show that (VYt(VX t ) _1 <7(i, X t )) t ^ T , is continuous. By assumption, a is a 

continuous function and it is well known that both processes (VX) _1 and X have continuous paths. 
VY is continuous because its dynamics is given as a sum of a stochastic integral of a predictable 
process against a Brownian motion (so a continuous martingale) and a Lebesgue integral with well 
behaved integrand. If two processes are versions of each other and one is continuous then they are 
in fact modifications of each other and hence Z has continuous paths. Now since Z has continuous 
paths, then the representation formula (I52p does not only hold dP <8> dt-almost surely but in fact 
holds for all t G [0, T] and P-almost all to G Q. Using that VY G Sq for some p > 2 (see Corollary 
13.11 and Proposition 13.21) . (VX)^ 1 , a(-, X) G Sq for any r > 2 and Holder's inequality, we conclude 
that Z G Sq for every q G [2,p). 

The property concerning the increments of Y is easy to prove since X, Y, Z G Sq for some p > 2. 
For < s < t < T, we have (recall that |/(-, ©(•))! < \f (■,&{■)) — /(-, 0, 0, 0) | + |/(-,0,0,0)j and 
that |/(-, 0,0,0)| is uniformly bounded) 

Y t -Y s = 0+ f f(u,@( U ))du- f Z u dW u , 

J s J s 

so using the assumptions and the Burkholder-Davis-Gundy inequality, we get for a generic constant 
C which may vary from line to line and some p > 2 

r r l 

E[\Y t -Y s \P] <CE f(u,e(u))du 

L is 

<c\t- s \^(i + ||xir s , + + \\z\f sE ) +e[( f \z u \^u) P ' 2 ] 

</ s 

<C\t-s\ p/2 . 

This in particular yields the applicability of Kolmogorov's continuity criterion to Y. □ 
The L 2 -regularity result 

We finish this section with the L 2 -regularity result for the control component Z of the solution of 
the time delayed FBSDE. Let it be a partition of the time interval [0, T] with N points and mesh 
size \ir\. We define a set of random variables via 



V 




p- 


+ 


y f z u m u 





Z% = E 



ZAs iTi 



, for all partition points ti, < i < N — 1. 
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The best square integrable JF^ -measurable approximation of t ^_ t . ft %+1 Z s ds is given by Z£., i.e. 



E 



U+i — t 



ti + l 



ZAs-Zl\ 



% Jt 



inf E 



ti+l — ti 



U+i 



ZAs - V\ 



(53) 



We associate the process (Zf) tG [o,T] to {Z£}i=o,--- ,n~i via ZJ = Z£ for t £ [ti,ti+\), Q < i < N — 1. 
Similarly, for the set of random variables {Z ti : ti G tt}, we associate the process (^f)tg[o,T] y i a 
= Z£ for t G [ij, tj+i), < i < JV — 1. The definition of the conditional expectation implies that 
for every i = 0, . . . , N — 1, we have 

E[|Z£j 2 ]-2E[Z£Z£]>-E[|Z£| 2 ], 

from which it follows that Z™ is the best ^-approximation of Z, leading to 



\Z-Z*\ 



K- 



< \\z -z*\ 



H 2 



0, as \ir\ 



0. 



Using Theorem 14.31 we are able to determine exp licitly the rate of conve r gence of the above limit. 
The following result extends Theorem 5.6 from llmkeller and Dos Reisl (|2010al ) to the setting of 
FBSDE with time delayed generators. 

Theorem 4.4 (L 2 -regularity). Assume that the conditions of Theorem \4 - 3\ hold for some p > 2 and 
assume further that a is ^-Holder continuous function in its time variable. Then 



max 

0<i<N- 



N-l r t i+1 

{ sup E[\Y t -Y t f}\ + Y,^\ * \Zs~Zl\ 2 d 

i ^u<t<t i+1 J ^5 l Ju 



< C\tt\. 



Proof. The result concerning the Y component follows immediately from Theorem 14.31 As for the 
result for Z, let us remark that since Z n is the best ^-approximation of Z over tt in the sense of 
(1531). it follows that 



Af-l 



iV-1 



N-l 



j2 e[ r +i \z s ~ zffds] <y, e [ I * +1 \ z * - z **i 2 H = e r +1 ^[\z s - z u \ 2 ]ds, 

where the last equality follows from the use of Fubini's theorem to switch the integration order 
(recall that Z G Sq for some p > 2). Theorem 14.31 allows to use (|52|) to rewrite the difference 
inside the expectation. We have Z s — Z ti = I\ + 12 + h with I\ = [VY S — VY ti ](VX ti )~ 1 a(ti, X ti ), 
h = VY S [{VX S )- 1 - {VX t .)- l }a{U,X ti ), h = VYsjVXs^lafaXs) - a{U X t .)} and s G [U,t i+1 ]. 

From the proof of part (ii) of Theorem 5.8 in llmkeller and Dos Eeisl (j201f)bh one obtains that 



N-l 



i=0 



ti+l 



|/ 2 | 2 ds + 



ti+l 



\h\ 2 ds 



< CM. 



The calculations that lead to the above result are quite easy to c arry out. They rely on known 
estimates for SDEs found for instance in Theorem 2.3 and 2.4 of llmkeller and Dos Reisl (|2010al ) 
combined with the fact that VY G S p for some p > 2. 
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To handle the term I± one needs to proceed with more care. Let us start with a simple trick: 



E 



\{VY s -VY u ){VX ti y l a(t u X u )t 



E 



^[\VY s -VY u \ 2 \T u ]\{\7X u y l o(t u X u )\ 2 ]. (54) 



Writing the BSDE for the difference VY S — VYf . for s & [ti, ij+i] we get for a generic constant C > 
that 



E 



|VY s -VY t .| 



< CE 

< CE 



U 



7T 



;(V/)(r,9(r)),(VG)(r)>dr| 2 + 



|(VG)(r)| 2 dr + 



7, 
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\VZ r \ z dr 

where we used the uniform boundedness of the derivatives of /, Jensen's inequality, Ito's isometry 
and proceeded to maximize over the time interval Combining the last line with (|54p and 

using the tower property, we obtain 



N-l 

E 

i=0 



rU+i 




/ E 


E 


hi 




N— 


i 



IVY, - VY fc 



T ti \(VX u y l a{U,X ti )f 



As 



< c M E 



i=0 





rU+i 


[(w 


hi 



ti+l 



(VQ)(r)| 2 dr + / \VZ r \ 2 dr)\(VX u )- 1 a(t u X ti 



N-l 



< |vr|E sup |(VX t )-V(t,X t )| 2 



0<t<T 



i=Q 



ti+l 

2, 



U 



|(VG)(r)| 2 dr + / |VZ r | 2 dr 



= |vr|E 

< CM 



sup |(VX t )-V(t,X t )| 2 (M / |(VG)(r)| 2 dr+ / |VZ r | 2 dr 

0<t<T v 



where in the last line we used the fact that VX, (VX) -1 , X G <Sq for every q > 2 and that VY, G 
T-£q for some p > 2 (in combination with Holder's inequality) to conclude the finiteness of the 
expectation. Combining this estimate with the ones for I2 and ^3 finishes the proof. □ 



Towards a time discretization of delay FBSDE 

Having established a path regularity result for FBSDE with time-delayed generators one can now 
start discussing a working numerical scheme. Given the nature of this class of BSDE, a time 
discretization would naturally require some decoupling technique to handle the backward-in-time 
feature of the equation and the backward-in-ti me feature of the delay . 



Applying the backward time discretization from iBouchard and Touzil (|2004l ) to (|33p -(|34p. we obtain 
for a partition tt : = to < h < . . . < tjy = T with step size Aj = ti + \ — ti 



Y tl = 9{X? N 
-W1 



Zl = E 



\ W U +1 - 
L A,: 



Wt, 



-Ytl +1 \Fti 



E 



Y t ^ +l \Tt] +A i f(t i ,&l] 



where 9£ = ( ^X^a^tj, t j+1 )), £ l£a y (fe, t j+1 )), £ Z^.a^tj, t j+1 )) 
j=0 j=0 3=0 
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This backward scheme cannot be implemented because in the computation of each Yp running 
backward from i = N — 1 to i = 0, we must evaluate Q w (ti) which depends on all Zf. running 
in forward direction j = 0, . . . , i. 

However, Bender and Denk ( 20081 ) propose for standard Lipschitz BSDEs a time discretization 
which mimics the Picard iteration technique for proving existence and uniqueness of BSDEs. Due to 
the fa ct that in each iteration step, one solves an explicit BSDE, the scheme from lBender and Denk 
(]2008l ) runs forward in time. The price to pay is to control apart from the error contribution 
of the time discre t izatio n the additional error arising from the Picard iterates (see Theorem 2 in 
Bender and Denk ( 20081 )). This idea adapts to equations (|33p -(|34 p by exploiting the fact that the 



solution (Y, Z) is obtained as a limit of (Y p , Z p ) as p goes infinity. Setting up (Y°,Z°) = (0,0) and 
then for p € No we have 

Y t p+1 = g(X T ) + ( T f(s,QP(s))ds- [ T ZV +1 dW s , t E [0,T] 
where Q p (t) = ( J° X t+V a x (dv), J° Y t p +V a z (dv), J° Z p t+V u z (dv)) . 

The discretization hereof is for p € No, initiated by setting (Y 71 "' , Z^^) = (0, 0), then iteratively for 
p > 1 and < i < N - 1 



YT 



7T,p+l 



z 



7T,p+l 



E 



E 



JV-l 

j=i 

w t .., - w t 



N-l 



A; 



j=i+i 



where 8^ = ( ^ ^ a x ([ tj , t j+1 )) , ^ Y t ^a y ([t jt t J+1 )) , ^ Z?/a z ([ tj , t j+1 ] 

j=0 j=0 j=0 

The proof of convergence for this time discretization scheme is left for future research. 
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